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Abstract. We present a brief overview of the regularizing transformations of the 
Kepler problem and we relate the Euler transformation with the symplectic structure 
of the phase space of the TV-body problem. We show that any particular solution of 
the TV-body problem where two bodies have rectilinear dynamics can be regularized by 
a linear symplectic transformation and the inclusion of the Euler transformation into 
the group of symplectic local diffeomorphisms over the phase space. As an application 
we regularize a particular configuration of the restricted circular TV+2 body problem. 
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Introduction 

In Celestial Mechanics, the iV-body problem has two types of singularities: collisions 
between two or more bodies, and escapes in bounded time. In order to study the 
behavior of the system close to singularities, it is a common procedure to transform it 
to another equivalent system that avoids the singularities by means of some methods 
called regularizations. There are a lot of regularizing transformations, but unfortunately 
it is not always possible to regularize an arbitrary singularity. For instance, infinite 
expansions at finite time produce esential singularities in the mathematical model that 
are not regularizable by topological or analytical methods known until now, and the 
same is true for some multiple collisions. 

Basically there exist two types of regularizations: analytic regularization formalized 
by Siegel and Moser [21], and surgery or topological regularization (also known as block 
regularization) discovered by Conley and Easton [8, 9]. In particular, it is well-known 
that rectilinear binary collisions, collisions between two infinitesimal bodies (in N+u 
problems) and triple collisions are impossible to regularize by the Easton method [19]. 
Marchal [19] has a very clear exposition about the classification of the singularities in 
the iV-body problem and their regularization (when it is possible). 

In this paper we deal with the analytical or Siegel's regularization [21] which is 
achieved by three ingredients: a local change of coordinates by means of some local 
diffeomorphism p : M — > M on the phase space, a scaling function g : M — >■ R 
that introduce a new fictitious time r by the relation ^ = g(w) and a set of initial 
conditions 0o = 0(0) of the flow which specifies the solutions that go to collision; since 
iV-body problems are Hamiltonian problems, the set of initial conditions is fulfilled by 
the conservation of the energy. It means that in Hamiltonian problems the analytical 
regularization process is performed on each fixed energy level H(x) = h. 

Thus, the process is as follows: 

• choose a fixed energy level H(x) = h and consider 

(H — h)(x) = 0, 

• apply the change of coordinates x = p(w) of the phase space 

[(H-h)op]{w) = 

• apply the scaling transformation 4^ = g(w) multiplying the last expression 

T=[g-(H-h)op}( w ) 

• the preimage of T(w, h) = generates the energy levels of the regularized system 
for each h G lmg(H) C R fixed. 

It is important to keep in mind that the aim of regularization theory is to transform 
singular differential equations into regular ones, controlling the velocity of the 
regularized system by the scaling time [7]. 
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For the one-dimensional Kepler motion with Hamiltonian function 

H(q,p) = ly 2 -^, xeW, y e T X R, (1) 
2 x 

it was already found by Euler that the introduction of a square-root coordinate u = y/x 
and a fictitious time r defined by dt = xdr reduces the Kepler equation of motion (1) 
to the equation of motion of a one-dimensional harmonic oscillator 

H = -v 2 + h u 2 , u e R*, v e T U R, (2) 

if h< [4]; where y = § , v = g and R* = (R - {0}). 

Generalizing this approach, Levi-Civita introduce its " .. .transformation du systeme 
qui donne lieu a des consequences remarquables..." in [17, 1907]. In his work Levi- 
Civita introduces a conformal transformation and exploits the symplectic structure of 
the complex plane (C, dz A dz) = (R 2 ,dy A dx). In fact, this regularization is made 
on the cotangent bundle T*(C*) where C* = G C : 2 ^ 0} is viewed as an open 
symplectic manifold. Levi-Civita regularization is achieved by the local diffeomorphism 
p:(T*(C*),u) (T*(C*),u) 



[z,w) 



,2 _w_ 

' ' 2|z| 



(3) 



and the time rescaling dt = \z\ 2 dr. The transformation above takes the Hamiltonian 
function 

H(q lP )= 1 - P 2 --^, P eT q C, (4) 

^ \1\ 



to 



fI =- w 2 + h\z\ 2 , zeC*, weT z C, (5) 



where p = w = and the symplectic form in the regularized phase space is 
u = dw A dz. The expression (3) is a contact transformation since it preserves the 
canonical Liouville 1-form a = wdz. If we denote the image of the local diffeomorphism 
by (q,p) = p(z,w) then 

pdq = wdz, (6) 

that is, p*(a) = a. As a consequence we have a symplectic (canonical) transformation; 
applying the exterior differential to both sides of (6) we obtain the symplecticity 
condition p*(oj) = u for the transformation. In 1913 Sundman introduced a 
transformation that maps the unitary circle in R 2 to the band — 1 < y < 1, and obviously 
this mapping does not preserves the area [10, pp 127-129]. 

Unfortunately, the procedure described above is difficult to generalize to the 3- 
dimensional case since the euclidean espace R 3 does not posses any complex structure. 
However, the Kustaanheimo-Stiefel's regularization [15], generalizes the Levi-Civita 
regularization to the four dimensional complex manifold T*C 2 (real dimension 8) and 
projects it onto some symplectic submanifold of real dimension six. In recent years, the 
K-S transformation using quaternions and the quaternionic algebra has gained much 
attention, from the works of Vivarelli [25], Volk [26], Vrbik [27], Waldvogel [28, 29], 
among others 
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1. Symplectic Structure of Regularizing Transformations 

In symplectic geometry, mechanical problems are represented by Hamiltonian systems 
(M,cu,X H ) on the phase space viewed as a symplectic manifold. The standard 
symplectic manifold is the cotangent bundle M = T*Q of the configuration space 
Q = (M n — A), and A is the set of the singularities of Xh and H. This manifold 
is provided with the canonical symplectic form uj = dp A <iq where q G Q and p G T* Q. 

In particular, problems on celestial mechanics are based on the Newtonian iV-body 
equations, Aiq = — ^ where 

N 

V(q) = G^^^, q=( qi ,...,q N ), ^1", 1 < i < N, (7) 
7T 'I' ~ 

and qi is the position of the i-th body, m ; its mass and M. = diag(mi/ n , . . . ,m,Nl n )- 
Depending on the value of n, we refer to this problem as the rectilinear or collinear 
problem if n = 1, the planar problem if n = 2 and the spatial problem if n = 3. In its 
Hamiltonian formulation the Hamiltonian function H : M — > M. is defined by 

if(q,p) = T(p)-V(q) (8) 

where T(p) = p T A^ _1 p is the kinetic energy. 

It is immediate that the set of singularities belongs to the potential function V(q). 
As we have said, it is not always possible to regularize any arbitrary singularity, however 
in this paper we are concerned with sigularities due to binary rectilinear collisions as 
the generalization of the rectilinear Kepler problem. To avoid this type of singularities 
we perform a regularizing transformation using a local diffeomorphism / : M — > M and 
a time rescaling g : M — > R. In some specific cases when it is desirable to preserve the 
fibers and sections of the cotangent bundle, the diffeomorphism and the time rescaling 
are applied to the base space / : Q — > Q and g : Q — > M. To obtain a symplectic 
diffeomorphism on M = T*Q it is used the properties of the cotangent lift of /. 

Definition 2 Let Q be an arbitrary differentiable manifold with cotangent bundle 
M = T*Q, and let f G Diff(Q) be any diffeomorphism over Q, we define the cotangent 
lift off by 

F:=T*f:M^M; F( Pl ) = p 2 , (9) 
where p { = (xi, &) G M, i = 1,2 and x { G Q, & G T*.Q, and (df Xl )*^ 2 = £i, 
Adittionaly we can see that 

(df Xl Y : T^Q ^ T^Q, 
so the restriction of FW* is the inverse mapping of (df x ,)*. 

Proposition 2.1 The cotangent lift F of any diffeomorphism f G Diff(Q) is a 

symplectomorphism, which means that 

F*uj = uj. 

where u is the canonical symplectic form on M. 
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The standard references where the reader can check the proof are [2, pp 487] and [1, pp 
180]. It is easy to show that the mapping 

T* : Diff(Q) Sp(M,w) 
/ ^F:=T*f. 

is a homomorphism of groups. 

In this way, it is possible to construct symplectomorphisms that preserve the 
structure of the cotangent bundle in the sense that they are fiberwise transformations. 
They form a subgroup of Sp(M) closely related to the set of generating functions on 
M. 

For iV-body problems in the plane it is a common procedure to identify the real 
plane with the complex numbers M 2 = C. Szebehely [23] has noted that in order to 
have a suitable regularizing transformation for binary collisions in the restricted plane 
3-body problem, the conditions 

z = f{w) and ^- = \f w ( w )\\ (10) 

must hold, where / : C — > C is a meromorphic function of the complex variable 
w = u + iv. 

The expression (10) is a fiberwise transformation which preserves the cotangent 
bundle as consequence of the cotangent lift of / : C — > C to T*C In such a case, 
the bilinear form uc = dp z A dz gives the symplectic stucture to T*C. Moreover, any 
fiberwise symplectic regularization of binary collisions in the N center problem has the 
form (10) [13]. This condition can be generalized for symplectic regularizations in higher 
dimensional spaces as it is exposed in [11] . 

As we have said, it is known by Euler that the transformation u = \fx and the time 
rescaling dt = xdr reduces the one-dimensional Kepler problem to the one-dimensional 
harmonic oscillator for h < 0. This transformation can be rewritten as x = u 2 with time 
rescaling dt = u 2 dr and it fulfills condition (10) when we restrict x > and / : M* — > M. 

In order to simplify calculations we add the coefficient 1/2 to the transformation 
and considering the cotangent lift we obtain 

U 2 V 2 

x = — , v = — , at = u dr. 11 

where y = ^§ and v = In what follows, we rename the variables x = q,u = Q,y = p 
and v = P to agree the standard notation of Hamiltonian mechanics. 

Definition 3 Let M = {x G 1R : x > 0} be the positive open ray and let V = T*J\f be 
its cotangent bundle. We define de Euler transformation £ : V — > V as the mapping 

(12) 

where Q e V and P G T£V. 

We restrict the domain of the Euler transformation to be an open manifold with 
boundary, in order to consider this transformation as a local diffemorphism. 
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Lemma 3.1 The Euler transformation £ : V — > V defined in (12) is a (local) 
symplectomorphism. 

Proof. We obtain the result in a straightforward way since the Jacobian matrix 



Definition 4 We call the Euler regularization of the collinear Kepler problem to the 
Euler transformation together with the rescaling function dt = Q 2 dr applied to the 
equation of movement of the Kepler problem. 

It is possible to consider the inclusion of the Euler transformation into the group 
Diffo(M) of local diffeomorphisms of any symplectic manifold (M, ui) containning a 
two-dimensional linear symplectic subspace V such that M = V © V^. 

We recall that a subspace V C E of some symplectic vector space (E, oj) of 
dimension 2n, is called symplectic if the restriction of the symplectic form u\y is injective 
(non degenerate). 

A well-known result about symplectic vector spaces that will be useful to understand 
the regularizing transformation applied to the circular iV-|-2 Sitnikov problem is the 
following. 

Lemma 4.1 Let (E, u) be a symplectic vector space and let V C E be a linear subspace. 
Then V is a symplectic subspace if and only if 



where V w is the orthogonal subspace to V with respect to the bilinear form uj. Moreover, 
V M is a symplectic subspace. 

The proof of this result is found in any book on symplectic geometry. Now, 
we procede to construct the regularizing transformation that we will apply to some 
symmetric N+2 body problems in the simpler cases: regularization of binary rectilinear 
collisions of the infinitesimals. 

Definition 5 The canonical inclusion of the Euler transformation into the group 
Diffo(M) of an open symplectic manifold with boundary (M,w) is the local 
diffeomorphism 




(13) 



is symplectic. 



□ 



E = V © V 



(14) 



% : V © V u -)• M, 



(15) 



such that 



(16) 
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We have the relation «o( = ^oj, therefore the following diagram comutes 

k 



M 



M 



where i(T*N) = V is such as in Definition 4. 

Lemma 5.1 The canonical inclusion is a local symplectomorphism G Spo(U,u) for 
U C M open subset. 

Proof. This fact is obtained straightforward from de direct sum V © V w , then the 
Jacobian matrix of the differential is 











^2(71-1) 



(17) 



where d£ G M.2x2 is the Jacobian matrix of the Euler transformation and htn-i) is th e 
identity matrix in 7W 2 (n-i)x2(n-i)- 

□ 



6. Some symmetric N+2 body problems 

Now, we must to characterize particular solutions of the N body problem where the 
Euler regularization is applied in a natural way. Since Euler regularization only considers 
the unidimensional (rectilinear) evolution of the colliding bodies we focus our attention 
to systems with N massive and 2 infinitesimal bodies and we call them N + 2 body 
problems. 

Definition 7 We say that a solution tp(t) = ip(ipo,t) of the spatial N body problem has 
a ^-symmetry around the line £cK 3 if there exists & G 50(3) such that the following 
fact holds: 

• for every t G (a, 0) and every sort S = (qi(t), • • • , <?jv(£)) the action of & on S is a 
cyclic permutation of order r > 1, 

• for every x G £ we have Mx = x. 

It is clear that the ^-symmetry applies to the whole phase space since this is valid in 
the configuration space for every t G (a, f3). This is equivalent to select a ^-symmetric 
initial condition ip(0) = tpo in the phase space M = T*(M? N — A) and follow the flow 

Proposition 7.1 Let tp(t) be a ^-symmetric solution of the spatial N body problem 
for t G (a, 0), around the fixed line £ctf and we consider the restricted N + 1 body 
problem attaching an infinitesimal bodies to the ^-symmetric solution. If the restricted 
body has initial conditions 

(#$) = (?)■ q> ° €R3 ' haa ' 0) (18) 
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such that q Q j G C and q® Ap® = 0, then the secondary evolves in rectilinear motion on 
the line C for t G (a, b) C (a, (3) . 

Proof. Since we are concerned with the evolution of the infinitesmal body with position 
g° G C, it is sufficient to show that piit) is parallel to C for t G (a,/3). 

By hypothesis (p (t) is a regular ^-symmetric solution of the N primary bodies 
around the line ££ for t G (a,/3). Without lost of generality, we can assume that the 
center of masses of the system is fixed at the origin and J?? = (0, 0, r), r G R is the 
vertical line in the 3 dimensional physical space. We supose also that the constant of 
universal gravity is G = 1. The ^"-symmetry implies that there exists a natural r > 1 
such that r\N and M r = Id, then & is a fixed matrix in 5*0(3) with components 

/ cos — sin ^ \ 



r r 



(19) 



-sin^ cos^ 

r r 

V l l 

Let s G N be the number of equivalent subsystems of the N body problem under the 
symmetry so N = rs. We can decompose the N body system in s partial subsystems 
with r bodies each one, in rearranging the subindices in the way 

{l,2,...,N)^{h, ■■■ ,l s ,2 1 , ••■ ,2 S , ■•• ,r 1; ■•• ,r s ) (20) 



such that 



and 



rh k := m kl = m k2 = . . . = m kr (21) 



q kj = Mq hj _ 1 =^ j ~ 1 q kl , (22) 
for k — 1, . . . , s and j = 2, ■ ■ • , r. Positions of each subsystem can be written as 

(9fci»?fc2,"--,?fcr) = (?fci»^9fci»-"«^ r ~ 1 ?fci)> fc = l,-"»s, (23) 
where g fc = (x k y k z k ). 



The Hamiltonian field for the infinitesimal body is 

mimi{qi - q { 



qj = —pi (24) 

771/ 



Rewritting equation (25) with reindexing (20) and expressions (21) and (22) we have 

P/=_m/ ^ P"^ | gf -^ gfcl |a ) (26) 

(27) 

We write qj = {xi,yj,zj) and pj = {p Xl ,Py I iPz I ) an d by hypothesis we have qj G C 
and qi Api = 0, it means that qj = (0, 0, zj) and pj = (0, 0,p 2j ). Finally, we note that 
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Y^-^Qki = (0,0,7*2^), we obtain the vector field as 

qj = (o, 0, — Pzi ) (28) 
\ mi J 



Pi 



0,0, -m/Vrmfc .^ \ . (29) 
which confirms that £ is invariant under the dynamics of the infinitesimal body. 

□ 

We assume known the solution </? : M x I — >■ M of the iV body problem, defined by 
ip(t) = tp((po] t) with ipo a ^-symmetric initial condition. Adding a second infinitesimal 
body with same conditions to those of Proposition (7.1), we obtain an iV + 2 body 
problem such that both infinitesimals have masses of the same order and they evolve in 
the vertical line for t e (a, b) C (a,/3). 

Without lost of generality, we can assume that the infinitesimal bodies have indices 
i = 1,2, with coordinates = (xi,yi,Zi) and pi = {j>x i -,P yi -,Pz^) and any element x e M 
is written as 

x = (xi 

) V\i Zl, ' " " ) l/N+2, z N+2iPx\iPy\iPz\i ' ' ' i Px N+2 > PyN+2^ Pz N +2) • 

Proposition 7.1 permit us to denote the position and momenta of each infinitesimal 
by (qi,Pi) = (0, 0, z iy 0, 0,p Zi ), i = 1,2, and their masses by m t and m 2 , the dynamics of 
both infinitesimals is given by the Hamiltonian vector field 

1 . „ zi - z k m x m 2 

zi = — p Zl , p Zl = -mi } rm k 

Til -i ^— » 



^ |^i -fe I 3 l^i-^l 



1 - ^2 — z k mim 2 

Z2 = — Pz 2 , Pz 2 = ~m 2 > rm k - - + 1 -, (30) 

m 2 ^ \z2-qk 1 \ ,i \zi-z 2 \ l 



with Hamiltonian function 

1 2 1 2 / v— rh k 



H(zi, z 2 , Pzi ,p Z2 ,t) = —p 2 n + — -p 2 Z2 - mi [Yl 



2m^ Zl 2m/ Z2 f^\7^i 1*1 -^'"Vl 

where q kj = M^^q^it), j = 1, ■ • • , r. Elements q^, for A; = 1, • • • , s, are representants 
of every cyclic subset under & which holds \zi — q kl \ = \zi — M^q kl \, for i = 1,2, and 
j = 1, • • • , r, as we have used in (30) to simplify the expression of the vector field. 

Remark 8 It is importante to note that the evolution of the N primary bodies is not a 
relative equilibrium in general; this is the case of the Hip-Hop and other more general 
solutions. As a consequence, is not always possible to reduce the dimension of the vector 
field. However, some configurations as the M -Circular Sitnikov problem [20] and the 
Sitnikov restricted N body problem [5] are reducible to one degree of freedom f. 

| In [20], Marchesin defines the M-Circular Sitnikov Problem where the configuration has N primary 
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Every solution of (31) will have a singularity due to collision when \z\ — z 2 \ — > at 
t — > b if b G (a, /3). We can regularize this type of singularities in order to extend 
the solutions of (31) for every t G (a, f3). Before to state the principal result of 
this section, we prove some technical lemmas which facilitates computations of the 
regularizing transformation. 

Lemma 8.1 The linear transformation Tb G Aut(M) with associated matrix B G 
Af6(Ar+2)x6(iv+2)(R), w/mc/j sends 

(z 1 ,z 2 ) ^ {z 1 -z 2 ,(l-fi)z 1 + iiz 2 ), fie (0,1/2] 

(Pz 1 ,Pz 2 ) H> (^P*! - (1 - ^Pz 2 :Pz 1 +Pz 2 ), 

and fixes all other components, is symplectic. 

Proof. It is sufficient to show that the reduced matrix 
/ 1 -1 \ 
(1-//) h 
fi -(l-fi) 

\ o oi i y 

is symplectic and a straightforward computation shows that B T JB = J. 

Since the masses of the infinitesimal bodies have the same order, it is possible to 
assign a linear relation in the form rri2 = c mi for some constant c G (0, 1]. Then the 
Hamiltonian function and vector field can be written in a more symmetric way. 



B 



(32) 



□ 



Lemma 8.2 The parameters m and e defined by 

mi + m<2 TTh\ — ni2 



m 



and 



mi + m 2 



and the time rescaling 
t = mt 

take Xh and H defined in (30) and (31) respectively to the form 



(33) 
(34) 



-Pzi, 



e) ^rm fc 

k=l 



-Pz 2 i P 



~2 



■d-«)E 



rrrik 



k=i 



Z\ — Zfa 

\zi - qk x \ 



m- 



+ m- 



\zi - z 2 \ 

1-e 2 
\z\ - z 2 \ 



(35) 



and 



H 



2(1 



\V Z1 



2 1 



EE 



k=l 




0=1 
1 - 



— m- 



Zi - z 2 



(36) 



bodies of mass m in circular relative equilibrium such that M = mN, and the infinitesimal in the 
conventional way. In [5], Bountis and Papadakis denotes the Sitnikov restricted N body problem the 
configuration with N — 1 primaries in circular relative equilibria. Here we follows the convention about 
N + v implies TV primary and v infinitesimal bodies as in [14]. 
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where H = H(z 1 , z 2 ,p Zl ,p Z2 ,i) , ' = j* and H = ±H. 

Proof. Direct substitution of the new parameters (33) into (30) and (31) 

Expressions (35) and (36) are obtained □ 

Theorem 8.3 The binary collisions between the secondary bodies of the symmetric 
N + 2 body problem with Hamiltonian function (31) are regularizables by the composition 
of a linear symplectic transformation A £ Sp(M) and the Euler regularization (i^, dt/dr) 
of the rectilinear binary collisions. 

Proof. First, we use Lemma 8.2 to work with normalized infinitesimal masses in such a 
way that Hamiltonian function (31) and vector field (30) are transformed to (36) and 
(35) respectively which depend on e and m as parameters. 

Let M be the phase space of the ^-symmetric N+2 body problem such that M is a 
cone m the total cotanget bundle M C T*M. 3 ( N+2 \ Since the evolution of the TV primaries 
is under a ^-symmetry, we can consider that C = {u £ M 3 | u = (0, 0, r), r £ R} is the 
symmetry axis of 

By Lemma (7.1), C is invariant under the evolution of the secondaries, then we 
are concerned with the third component of their coordinates q^ = (0,0, z/J, and 
pj i = (0,0, p Z[ ), for i = 1,2. Since the indexing of coordinates will be tediuous we will 
assume that I\ = 1 and I 2 = 2 and the primaries will have coordinates qi = (xi,yi,Zi) 
and pi = (Px^Pyi^Zi), for 3 < i < N + 2, which permits to express a single point x £ M 
in the form 

x=(x 1 ,y 1 ,z 1 ,X2,y2,Z2,x 3 ,---,z m ,p xl ,---,p Zm ), m = N + 2. (37) 

We select M to be the cone which holds z\ > z 2 such that the infinitesimal masses are 
in relation mi > m 2 . 

Consider the transformation Te £ Aut(M) which permutes the coordinates with 
indices 

1 2 3 4 5 6 7 ••• i ■■■ n + 3 n + 4 n + 5 n + 6 
5 6 1 7 8 3 9 •■■ i + 2 •■■ 2 n + 5 n + 6 4 

where n = dim(M)/2, in such a way that the determinant of the associated matrix E 
is unity. 

Matrix E have the following properties 

• E T = E~\ 

• E T JE = J. 

Define A = EB where B is the associated matrix of transformation from Lemma 8.1 
and A £ Sp(M) since E,B £ Sp(M). Then A sends any element x £ M to 

A ■ x = (zi - z 2 ,/J>p zl - (1 - n)p Z2 ,x il} ■ ■ ■ ,p ij ) T , /IE (0,1/2], (39) 

where [i = ^ and ■ ■ • , ij) is an even permutation of the other 6(iV + 2) — 2 indices. 

First two components of (39) define a symplectic subspace (V,w|y) of the phase 
space (M,uj). The other elements, define another symplectic subspace (V", uj\yu) which 
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is w-orthogonal such that M = V © V w . Singularities due to binary collision of the 
secondary bodies belong to subspace (V, o>|y) and this happen when first component 
goes to zero. 

Now, we can apply the canonical inclusion of the Euler transformation to A in the 

form 

Q = (i £ Oi4)x, Q,x g M, 
in such a way that the transformation 

p:=(i^oA)- 1 (40) 

and the rescaling time dt = 2/i(l — p)Q\dr regularizes the binary collisions of the 
system (35) and (36) by the analytic method. Since the set Sp(M, u) is a group under 
composition, we have immediately that transformation (51) is also symplectic on the 
manifold (M, u) and the following diagram commutes 



M M M 




The regularized phase space will have the form 

%:V©V w ^ft©V w . (41) 

where 1Z is the symplectic subspace V under the Euler regularization. 

Denoting by Q = (Qi, Q 2 , - • • , Q n , Pi, P2, • • • , P n ) G M we get in local coordinates 
the expression for x = p(Q) given by 

Zl = Q 2 + p 2l =/iP2 + ^, (42) 

Z2 = Q2~^, p Z2 = (l-fi)P 2 -^-. (43) 

All other components obey the rules of indexing of x given in (37) and permutation 
(38). Substituting (42) and (43) into (35) and (36) and applying the time rescaling with 
fi = and (1 — fi) = we obtain the Hamiltonian function in the form 

r = \ (^-PlQl + *?) " l ~^Ql WQuQ*) +h]-(l- 6 2 ) 2 m. (44) 
where 

v(n) ( (j+j)™* (l_-e)m^_\ 

' ttf^\\Q2+ l -^Qi-Q kj \ Q,-^Q\-Q k . 3 \)' 
Q = ((Qi,Q 2 ,o,o,o,o),(Q ll ,---,g Sr )), 
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and Qk 3 — Qkj, k = 1, • • • , s, j = 1, • • • , r, are the positions of the N primaries. Since p 
is locally symplectic on the open manifold with boundary (M, u) the new Hamiltonian 
vector field can be obtained directly from the regularized Hamiltonian function 

r:Ixl4l (45) 
(Q,h) ^r{Q;h) (46) 

which depends on the parameter h 6 R. The regularized Hamiltonian vector field X^ h 
which also depends on the energy level h has the form 



Q'i 


= Pl^ 




Q'2 


— 4 Q1P2, 




P[ 




'pi-4[V 1 {Q) + h]-2Q 1 


P 2 


-^Ql- 


8Q 2 - 



(47) 



Expressions (44) and (47) are free of singularities due to collisions between the 
secondary bodies. In this way, we are extended the Hamiltonian system to the set 
A := {zi = z 2 } which is a subset of the boundary dM. 

□ 

Remark 9 The energy levels H(x) = h are mapped to the zero set of Y and we will 
denote them by 

S ft = {Q G M r(Q; h) = 0, he lmg(H) C M} (48) 

regardless if this is the energy level in the original system or in the regularized one. 

The Hamiltonian vector field Xr h is valid only on the energy level for every 
h G Img(if) fixed. 

Examples of this type of systems are the circular collinear N + 2 and 2N + 2 
problems showed in Figure 1. Other examples are constructed with 2N massive bodies 
in a Hip-Hop solution and 2 infinitesimals bodies on the line determined by the angular 
moment of the system as the reader can see in Figure 2. 

Now, we proceed to study a special case of the iV+2 body problem with many 
symmetries. We called this problem the circular iV+2 Sitnikov problem [12] since this is 
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^ 










^ 








l — 7 



Figure 2. HipHop-collinear 2 AT + 2 problem. 



a generalization of the circular N Sitnikov problem [20, 5] which contains an additional 
infinitesimal body. 

10. The Circular N+2 Sitnikov Problem 

For an application of the symplectic regularization, we select a special configuration of 
the restricted N + 2 body problem. This particular configuration has N massive bodies 
with masses 777.3 = ■ ■ ■ = mN +2 — in relative equilibrium evolving in circular orbits 
on the vertices of a regular iV-gon around their center of masses. The system has two 
infinitesimal bodies that evolve on the perpendicular straight line which passes across the 
center of masses of the massive bodies. The massive bodies are called primaries and the 
infinitesimal bodies are known as secondaries. The problem consists in determining 
the evolution of the secondaries under the attraction of primaries with Newtonian 
gravitational potential. In general, the secondaries have different infinitesimal masses 
m i 7^ m 2 and without lost of generality we can assume that m 2 < m\ <C -h. 
Let Q be the configuration space defined by 

Q = {(91,92) eMV ><? 2 }, 

where g« is the position of the body with mass 777j for % = 1, 2. The potential function 
V : Q — > M of the circular N + 2 Sitnikov problem is 

V(qi, q 2 ) = , s + -, s + , 

V 9i + r V ?2 + r ?i _ 92 

and the Hamiltonian function H : T*Q — > K will be 

^(q,p) = ^p T -M~ 1 p-^ / (q), (49) 

where q = (91,92) is the vector of positions, p = (pi,p 2 ) is the vector of conjugate 
momenta and M. = diag(mi,m2) is the matrix of masses. The constant of universal 
gravitation is G = 1 and r is the radius of the circunference which contains the vertices 
of the iV-gon and must fulfill the condition [18] 

v 1 

2w 2 r 3 = Gm V— — — , N = 2u + 1, ,iseN 
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7= 

where w is the angular velocity and m is the mass of each primary in the circular relative 
equilibrium. Since G — 1 and m — vp, and considering w = lwe obtain the relation 

o 1 \ 1 o 

r = > or r 

2N 2- sin (^) 

whenever i/ is odd or even respectively. 

Equivalent expressions are found by Bountis and Papadakis in [5] in the N+l 
Sitnikov problem where the value for r is given by 

r = \csc 

and the masses of the primaries are m = with 

r- t A sin 2 9 cos(^ + 1 - i)9 

*r~i sin [v + 1 — i) 9 

i=2 v ' 

Marchesin [20] in contrast, fix the radius r = | and by a suitable rescaling studies the 
effect of the variation of primary masses on the period function T(h) = T(h;m). In 
general, a suitable change on the angular velocity and the masses of the primaries let 
us to normalize the radius to r = 1 (which is not the case in this paper). 

Applying Lemma 8.2 we will write the masses of the secondary bodies as mi = 
m(l + e) and rri2 = m(l — e) and the rescaling time t h- >■ mt will produce the reduced 
masses a = 1 + e and — 1 — e. Denoting by r = the radius of the circunference 
for the iV+2 Sitnikov problem, the potential function now depends on the number of 
primary bodies as a parameter, then K: QxN->l becomes 

1 + e 1 - e 1 - e 2 

V(q u q 2 ] N) = - + - + m- 



and the Hamiltonian function becomes 

1 9 1 9 1 + e 1 — e 1— e 2 , . 
H(q,p:N) = — -pi + — -pi , - m . (50) 

It is importante to note that the angular velocity of the primary bodies is no more 
unity w ^ 1 due to the time rescaling t — » mt, however this fact is not relevant in the 
rectilinear problem. 

Remark 11 The symmetry (qi,Pi,q2,P2,t) l— > (q2,P2, ft; Pi, ~ e ) restricts the analysis 
to non negative values of the parameter e. 

Let M = T*Q be the phase space of the Hamiltonian system TL = (M,uj,Xh) 
associated to the problem, where u = Y2i dpi A dqi is the standard symplectic form on 
M. A = {(qi,q2,Pi,P 2 ) £ ^ 4 |ft — ft} is the set of singularities of H(q, p; N) due to 
collisions and it is easy to see that A = dM. 
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Figure 3. The circular N + 2 Sitnikov problem for N = 3 and N = 8. 



The Hamiltonian vector field Xjj in local coordinates is as follows 



1 



qi 



Q2 



1 + e 
1 

1 - e 



Pi, 



P2, 



Pi 



P2 



(qf + r 



e)qi 1 
■ — m 



NJ" 



+ m 



NJ 



(?1 - ?2)' 



The evolution of both secondaries is restricted to the perpendicular line that passes 
by the center of masses of the primaries. The symmetries of the problem keeps the 
secondaries on the perpendicular line and since the angular moment is null there is not 
scattering at collisions. 



11.1. Regularization 

To avoid the singularity in both, the Hamiltonian function and the vector field Xh, we 
perform a symplectic regularization. In order to extend analytically the equations to 
the hyperplane q\ = q2 we apply the transformation p : M — > M defined by 



Qi = Q2 + 

and the time rescaling 
dt 1 1 



4 

1 + e 



4 



Ql 



Pi 
Pi 



+ e p + Pl 

—P2-^ 

2 2 Qi' 



-Ql 



IS 

r 



If we write z = (Qi, Q 2 , Pi, Pi) and \i 



( M (l - »)PiQ{ + Pi) - 16^(1 - fif 



(51) 



(52) 



the regularized Hamiltonian function 



m 



-2/i(l - fi)Ql 



4(1 -fi) 



4/i 



(2Q 2 + fiQlY + Ar 2 N 



(2Q 2 -(l-fi)QlY + Ar% 



the 



We denote T h (z,fi;e) = T(z, p; e, h), and we call to the triplet (M,u,X r 
regularized system, where M = T*QU{q\ = q 2 } and X-p h is the regularized Hamiltonian 
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field 



Q 



dP 



9l\ 
dQ 



In local coordinates we get 

Q[ = Pi, 



P-2 



li(l-fi)QfP 2 



H(l - /i)Q! (p| - 4 [V(Q U Q 2 ) +h}- 2Q 1 ■ ^) 



-2Mi-/i)g?-g. 



Computing the partial derivatives we obtain 

2Q 2 -(1-^)Q 2 



^ = 8/i(l-/i)Qi 



dV 



[(2Q 2 -(l- M )Q2)2+44]2 [(2Q 2 + M Q?) 2 +4r2 r ]3 
M(2Q 2 -(l-M)Qf) , (1-m)(2Q 2 +/^) 



[(2Q 2 -(1-m)Q?) 2 +4t^]2 [(2Q 2 +MQ?) 2 +4r2,]7 

and arranging equivalent terms in the expression 



IV + 2Q X ■ 



dV 
dQi 



16 



M (4Ql-2(l- M )Q 2 Q 2 +4r 2 r ) 



[(2Q 2 



(l^)Q 2 ) 2 +4r2,]r + 



(l-/x)(4Q 2 +2 M Q 2 Q 2 +4r 2 f ) 
[(2Q 2+M Q 2 ) 2 +4r 2 ,]2 



we obtain the vector field as 
Pi, 

Ki-m)Q?p 2 , 
-/i(l - 



Qi 
Q' 2 

P[ 



p2 
r 2 



Ah - 16 



(l-/x)(2(2Q 2 +MQ 2 )Q2+4r 2 f ) 



[(2Q 2 +MQ 2 ) 2 +4r iv 



2 if" 



M2(2Q 2 -(l-M)Q 2 )Q2+4r| r ) 
2 -(l- M )Q 2 ) 2 +4r 2 ,]7 



P 
r 2 



-16/i(l 



[(20 



M2Q 2 -(i-A')Qi) 



+ 



(l-p.)(2Q 2 +fiQ\) 



(53) 



(54) 



(55) 



(56) 



[(2Q 2 -(l-M)Q?) 2 +4r 2 r ] 3 [(2Q 2 +/,Q 2 ) 2 +4r 2 r ] 

Although the form of the new Hamiltonian function and the vector field are 
complicated, the advantage is that they are regular in M := M U A. 



quite 



11.2. Symmetries. 

The regularized Hamiltonian function has a symmetry in Pi and P 2 that reflects the 
symmetry with respect to the fictitius time r in the way 

(Qi, Q 2 , Pi, P 2 , r) ^ (Qi, Q 2 , -Pi, -P 2 , -r). (57) 

It is a generic property of mechanical systems. The symmetry in the Q\ variable 
is fictitius due to the transformation Qf/2 = q\ — q 2 . Finally, applying the change 
Q 2 H- — Q 2 it changes the values of (1 + e) H- (1 — e) and viceversa. 

Theorem 11.1 The regularized Hamiltonian system (M, u,X-p h ) is symmetric with 
respect to the hyperplane Q 2 = if e = 0. Moreover, if e = 0, the symplectic plane 

Si = {(Qi, Q 2 , Pi, P 2 ) e M\Q 2 = P 2 = 0} 

is invariant under the flow of the regularized Hamiltonian vector field X-p h ■ 
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Proof. Using the Hamiltonian function I\ and substituting Q 2 — > ~Q 2 it remains 
invariant if 



The last identity has as trivial solution — and this holds if and only if e = 0. 

In order to prove that S\ is an invariant plane under the flow we consider Q 2 = 
for every r G I C R. By hypotesis e = and consequently fi = |, then the fourth 
equation in (56) implies P 2 = and therefore P 2 = constant. Additionally, Q 2 = 0, but 
we know that /i(l — fi) 7^ and Qi is not identically zero. Then P 2 = and we have 
the reduced system 




Qi = Pi, P[ = -8Q1 

Q' 2 = 0, P' 2 = 0, 

where a = Ar^- Consequently, S\ is an invariant plane under the flow 0(t) of the 
Hamiltonian vector field XY h . 

□ 

It is known that Hamiltonian systems (M, u,Xh) which have invariant symmetry 
planes can be reduced to systems restricted to the invariant plane. In fact, each invariant 
plane corresponds to some symplectic subspace and vector fields restricted to symplectic 
subspaces can be locally integrable. In this example, the flow <Ph{ t ) of the Hamiltonian 
system restricted to the symplectic subspace Si is equivalent to have the secondaries' 
relative barycenter miqi + m 2 q 2 = at the origin. 

Definition 12 We define the symmetric circular N+2 Sitnikov problem to the 

Hamiltonian system (M, u, Xh) where e = and the initial conditions are symmetric 

It means that po = pi(0) = —^2(0) and qo = qi(0) = —q 2 (0). we have the following 

Corollary 12.1 The symmetric circular N +2 Sitnikov problem for m ~ 0, is integrable. 

Proof. It is an immediate consequence of Theorem 11.1. Since the initial conditions 
are gi(to) = -^2(^0) and pi(t ) = -p 2 (t ) and e = then Q 2 (r ) = and P 2 (r ) = 0. 
Additionally, Proposition 11.1 implies that Si is an invariant symplectic plane then 
Q 2 (t) = and P 2 {r) = for all r G / C M where / is its domain of definition. 

Therefore, the symmetric circular N+2 Sitnikov problem is a Hamiltonian system 
with one degree of freedom. It has as regularized system to (M, u, Af ) with M = T*~R* 
and regularized Hamiltonian function 

f = -Pl-AQ\ ( - 1 +-1 -m. (59) 

where a = ^r^. This is a first integral for the reduced Hamiltonian system when T = 

□ 

The vector field Ap in local coordinates is as the two first lines in (58) and the level 
curves are show in Figure 4 
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Proposition 12.2 The symmetric circular N+2 Sitnikov problem has the following 
dynamics: 

• If h < the solutions are periodic orbits where the secondary bodies collide at the 
origin of coordinates. 

• If h = the system has a parabolic solution with the escape of both secondaries with 
null velocity when they reach the infinity. 

• If h > the solutions are hyperbolic orbits with escape of both secondaries in oposed 
directions and with positive velocity at infinity. 

Proof. We verify this fact directly from the Hamiltonian function of the original 
system. Substituting e = 0, qi(t ) = — 52(^0) an d Pi(t ) = — ^2(^0) hi (50). Defining 
p := pi = —p2 and q :— qi = —q2 and fixing H = h we obtain 

- = - 2 - 1 _ ™ egg) 
2 2 P ^+r% 4g- 1 ' 

The maximum distance from the origin that secondaries can reach is for p = and then 
h 1 m 



2 vV+^i 4 <? 

has a finite real solution q > for every fixed h < 0. It means that the evolution 
is bounded and extending the solutions beyond collisions with the regularization, the 
solutions are periodic orbits with elastic bouncing at collisions. 

On the other hand, if h > this approach does not apply. For this case, we solve 
(60) for p = q to obtain 



q = ±Jh+-=== + -, (61) 




and we obtain the escape velocity by the limit 
lim q = ±\//z. 

q— ¥00 

Since we deal with the symmetric problem, we are concerned only with positive values 
for q, p and h. Negative values are associated with the other secondary body. 
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For h = the limit lim^oo q = implies that the bodies escape to infinity with 
cero velocity, which confirms the parabolic orbit. 

Finally, for h > we have lim^oo q > and the solutions are hyperbolic orbits, 
where the secondary bodies escape to infinity with positive velocity. 

The dynamics is as follows: secondary bodies start its evolution at infinity from 
oposed sides of the plane where the primary bodies evolve. Secondaries approach to 
massive system symmetrically to collide at the origin with ellastic bouncing and escape 
to infinity in oposed directions. 

□ 
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